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ABSTRACT 

Context. The mass-metallicity relationship (MMR) of star-forming galaxies is well-established, however there is still some disagreement with 
respect to its exact shape and its possible dependence on other observables. 

Aims. We measure the MMR in the Galaxy And Mass Assembly (GAMA) survey. We compare our measured MMR to that measured in the Sloan 
Digital Sky Survey (SDSS) and study the dependence of the MMR on various selection criteria to identify potential causes for disparities seen in 
the literature. 

Methods. We use strong emission line ratio diagnostics to derive oxygen abundances. We then apply a range of selection criteria for the minimum 
signal-to-noise in various emission lines, as well as the apparent and absolute magnitude to study variations in the inferred MMR. 
Results. The shape and position of the MMR can differ significantly depending on the metallicity calibration and selection used. After selecting 
a robust metallicity calibration amongst those tested, we find that the mass-metallicity relation for redshifts 0.061 < z < 0.35 in GAMA is in 
reasonable agreement with that found in the SDSS despite the difference in the luminosity range probed. 

Conclusions. In view of the significant variations of the MMR brought about by reasonable changes in the sample selection criteria and method, 
we recommend that care be taken when comparing the MMR from different surveys and studies directly. We also conclude that there could be a 
modest level of evolution over 0.06 < z < 0.35 within the GAMA sample. 
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1. Introduction 

The mass-metallicity relationship (MMR) describing a correla- 
tion between stellar mass and gas-phase m etallicity in galaxies 
was first reported by Lequeu xet alj d!979l) . Looking at spectra 
of H ii regions for a small sample of irregular and blue compact 
galaxies, they found that metallicity and stellar mass were cor- 
related over two orders of magnitudes in mass. Subsequently, 
with the advent of large spectroscopic surveys of galaxies, the 
reality of the MMR has been co nfirmed and establish ed many 
times since, for both stellar (e.g., Gallazzi et al. 2005; Mendel 
2009I) and gas-phase (e.g JTremonti et alJ2004l : IZahid et alJ 



etal. 



|201 1 ) metallicities. 

Several mechanisms have been proposed to explain the 
origins of the MMR. These include outflows, with enriched 
gas being retained preferentially by high-mass g alaxies due to 
the increased depth of their potential well (e.g., lLarson|[l974t 
iKobavashi et al J 12007; Spit oni et alj|2010t) . Another possible 
mechanism is the interplay between the outflow of enriched gas 
and infall of pristine g as from the inter-galactic medium (e.g., 
iFinlator & Dav e 2008). The downsizing phenomenon, wherein 
stars found in higher-mass galaxies today formed rapidly at ear- 
lier epochs suggests the surrounding gas should be enriched 
early. Stars found in low mass, low redshift galaxies, however, 



have lower cr-element abundance r atios, suggesting th ey formed 
slowly over longer periods (e.g.. ICowie et al] 19961 Kodama 



et al. 



2004tlAbraham et alj|2005t iThomas et alj|201ob . Indeed 



iGarnetj ( 12002 ) suggest that the MMR can be explained by the 



fact that low mass galaxies have higher gas fractions than high 
mass ones , implying that they are still converting gas into stars. 
Moreover. iRodrigues et alj d2012h show that galaxies at redshift 
Z ~ 2.2 have higher gas fractions than their local counterparts, 
suggesting that evolution of the MMR should be expected un- 
der this scenario. A fourth possible mechanism is a stellar mass 
dependent or star formation rate (SFR) dependent evolving ini- 
tial mass function ( IMF), which has also been proposed as a 
viable scenario e.g., Koppen et al. 120071 : IWilkins et al. 2008a; 
I Wilkins et~aT]l2od8bl:ISpitoni et alJl2010l:TGunawardhana et all 
1201 lHFerreras et alj|2012h . This would also impact the rate and 
level of metallicity enrichment with galaxy stellar mass. 



* cfoster@eso.org 



The Sloan Digital Sky Survey (SPS S) has allowed for the 
unequivocal confirmation of the MMR (Tremo riti et al.| [2004. 
hereafter T04). The method used in lT04l is based on simultane- 
ous fitting to all available emission lines using stellar population 
models. While using all the available lines to make a best esti- 
mate, this approach is more sensitive to noise. This method is ul- 
timately tied to the reliability of the models used, and artificially 
produces quantised output metallicities as finite values of metal- 
licities are modeled. lT04l provide a conversion from the popular 
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R23 = ([On]/l3727 + [Orn]/U4959,5007)/I^8 strong emission 
line ratio to their calibr a tion. 

Kewley & Ellison (I2008L hereafter KE08) have used the 
SDSS data to compare the inferred metallicities using various 
calibrations. Empirical conversions between the different metal- 
licity calibrations are provided. They find that the absolute posi- 
tion and shape of the MMR can vary significantly from one cali- 
bration to another, emphasising the need to use similar metal- 
licity calibrations when com paring different samples. In deed, 
KE08 and other authors (e.g. lLopez-Sanchez & Esteban 2010; 
Bresolin et al. 2009; Moustakas et al.l 120100 have shown that 
metallicities d erived using calibrations based o n photoionization 
models (e.g., lMcGaughlll99Tt IKD02I : IKK04 tend to be sys- 
tematically 0.2-0.4 dex higher than metallicities derived using 
the direct method (i.e., using a direct estimat ion of the electron 
temperature) or calibrations based on it (e.g., PP04; Pilyugin & 
Thuan T2005|). In addition, depending on the metallicitv calibra- 
tion used, KE08I found that the turnover of the MMR depends 
on the aperture covering fraction. To avoid such effects, they 
recommend a minimum redshift of z = 0.04 for the SD SS 3" 
aperture fibre spectra. Similarly, iMoustakas et al.l (1201 ll) find a 
measurable dependence of the MMR on the metallicity calibra- 
tion used and aperture biases using a sample of ~ 3000 galax- 
ies from the AGN and Galaxy Evolution Survey (AGES). They 
also find that the MMR is affected by contamination from Active 
Galactic Nuclei (AGN). Hence, differences induced by metallic- 
ity calibrations, aperture effects and AGN contamination need 
to be taken into account when comparing results from different 
surveys and studies. 

Several groups have studied the redshift evolution of the 
MMR and sometime s find conflicting results , particularly at 
lower redshifts. Bo th lLara-Lopez et alJ d20 09a.b) and Savaglio 
et al. d2005h detect evolution of the MMR up to z = 0.4 using 
SDSS, and over 0.4 < z < 1.0 from a combined sample from 
the Gemini Deep Deep Survey (GDDS) and the Canada-France 
Redshift Survev (CFRS), respectively. In sharp contrast, Carollo 
& Lilly (|2001D do not find significant evolution when comparing 
the MMR in a sample of galaxies at 0.5 < z < 0.7 from CFRS 
to that locally observed. Yet, several groups independently de- 
tect signi ficant evolution at comparable redshifts 0.4 < z < 1.0 
(see e.g.. iKobulnicky & Kewlevl 120041: iMouhcine et all [2006: 
iRodrigues et al ll2008l:lMorelli et alJfcom Pushing the redshift 
barrier further. lErb et al.l d2006l) analyse a sample of 87 UV se- 
lected star forming galaxies to show that galaxies at high red- 
shift (z ~ 2) not only have lower oxygen abundances, but also 
higher gas fraction (~ 50 per cent). At redshifts z > 3, Maiolino 
et al. (120081) select a variety of independent ratios and calibra- 
tions to obtain metallicities and detect an evolution of the MMR 
and its slope in a sample of 9 star forming galaxies observed 
with ESO-VLT. They claim that no hierarchical formation simu- 
lation reproduces the MMR at z ~ 3, w hile monolithic col lapse 
simulations do. At even higher redshift. [Laskar et al.l d201 ll) ob- 
tain a sample of 20 gamma ray bursts at redshifts 3 < z < 5 and 
demonstrate that the MMR at high redshift is significantly offset 
to lower metallicities than that at z < 3. 

Another fundamental observable parameter that has proven 
to be important in order to disentangle which scenario(s) is (are) 
responsible for the MMR is the SFR. Here again, there are some 
discrepancies i n the various fi ndings about the dependence of the 
MMR on SFR. lEllison et all d2008l) select galaxies in the SDSS 
to study the dependence of the MMR on the specific SFR (i.e., 
SFR per unit mass, hereafter SSFR). They observe that galax- 
ies with higher SSFR for their stellar mass have lower metal- 
licities (log(0/H) is as much as 0.2 dex lower). More recently, 



iMannucci et al.l d2010l hereafter M10) have also used SDSS 
galaxies to study the inter-dependence between mass, metallic- 
ity and SFR. They claim that there exists a 'fundamental rela- 
tionship' or three-dimensional curved surface in stellar mass, 
SFR and metallicity space where all star forming galaxies lie. 
They also find that higher SFR galaxies (and lower metallici- 
ties) are selected at higher redshift, leading to an apparent evo- 
lution of the MMR . lMlO l demonst rate that high redshift g alax- 
ies from lErb etal] d2006l) . but not IMannucci et alJ d2009l) fol- 
low this 'fu ndamental relationship' . Sim ilarly, lLara-Lopez et al. 
d2010al) and lLara-Lopez et all d2010bl hereafter LlOb) discov- 
ered a relationship which is best described as a "fundamental 
plane" b etween SFR, metallicit y and stellar mass (r ecently up- 
dated bv lLara-Lopez et al. 2012). In contrast to lMlOL they do not 
find any curvature in three-dimensions or any evolution of the 
fundamental plane up to z ~ 3.5, and detect a shallow but posi- 
tive correlation between SFR and metallicity such that high SFR 
galaxies tend to have higher metallicities. In general, they ob- 
serve that metallicities are lower, SFRs are higher and the mor- 
phology of galaxies shows a higher fraction of late-type galaxies 
at hig h redshift comp ared to locally. Using semi-analytic mod- 
els, ^tes_etal] (j20T2) reproduce the 3-dimensional surface the- 
oretically. They explain that low SFR, high mass galaxies have 
typically exhausted their gas reservoirs in a recent major merger, 
hence preventing further star formation. Subsequent inflow of 
metal-poor gas then dilutes the gas around these galaxies with- 
out significant further star formation. 

The focus of this paper is to present the MMR using data 
from the Galaxy A nd Mass Assembly survey (GAMA, see 
iDriver et al.ll20liTPl and study its dependence on sample se- 
lection. A companion paper (Lara-Lopez et al., in preparation) 
presents the joint distribution of GAMA and SDSS galaxies 
in the three-dimensional SFR/stellar mass/metallicity parame- 
ter space. GAMA is a good complement to the previous SDSS- 
based studies as it probes to fainter magnitudes than SDSS, in 
principle allowing for a range of lower stellar masses to be stud- 
ied at a given redshift, and to higher redshift at a given stellar 
mass. The paper is structured as follows. The data and sample 
selection are described in Section [2] Sections [3] and [4] present 
our analysis and highlight our results, respectively. A summary 
and our conclusions are given in Section|5] 

Throughout we assume a cosmology given by Q.m = 0.3, 
Q. A = 0.7 and H = 70kms-' Mpc" 1 . 



2. Data 

2.1. GAMA survey 

Upon completion, the GAMA survey (see IDriver et al. l l20Tll) 
will include data from many telescopes (AAT, VST, VISTA, 
Herschel, ASKAP, and GALEX) to create a database ranging 
from radio to ultraviolet wavelengths for roughly 400,000 galax- 
ies over 360 square degrees. Optical spectra from AAOmega on 
the AAT have been obtained for ~ 150,000 galaxies over three 
regions covering 144 square degrees in total as part of GAMA 
phase I. The GAMA Phase-I spectroscopic completeness is 98% 
to an r-band Petrosian magnitude of 19.4 over two of the regions 
and r ~ 19.8 over the third region. GAMA Phase-II will include 
a larger volume and will be complete to an r-band Petrosian 
magnitude of 19.8. Because of its depth, it is a good complement 
to SDSS probing to lower stellar masses and higher redshift. This 
work uses GAMA Phase-I data. 



http ://w w w. gama- survey, org/ 
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-1.5 -1.0 -0.5 0.0 0.5 

log( [NII]X6584/ Ha ) 



Sample 


Miiin 








(1) 


(2) 


(3) 


(4) 


(5) 


1 


0.070 


0.085 


-18.58 


-19.73 


2 


0.085 


0.100 


-18.97 


-20.18 


3 


0.100 


0.115 


-19.30 


-20.57 


4 


0.115 


0.130 


-19.60 


-20.90 


5 


0.130 


0.145 


-19.86 


-21.20 


6 


0.145 


0.160 


-20.10 


-21.46 


7 


0.160 


0.175 


-20.33 


-21.70 


8 


0.175 


0.190 


-20.53 


-21.93 


9 


0.190 


0.205 


-20.72 


-22.13 


10 


0.205 


0.235 


-21.07 


-22.32 


11 


0.235 


0.265 


-21.39 


-22.67 


12 


0.265 


0.295 


-21.67 


-22.99 


13 


0.295 


0.325 


-21.93 


-23.27 


14 


0.325 


0.350 


-22.13 


-23.53 



Table 1 : Definition of the selected volume limited samples (col- 
umn 1). Redshift limits are given in columns 2 and 3, while ab- 
solute Petrosian magnitude limits can be found in columns 4 and 
5 for each sample. 



Fig. 1: Selecti on of star-forming galax ies (red) using the pre- 
scripti on from iKauffmann et alJ (I2003L dashed line) BPT dia- 
gram (Baldwin et alJll9811) . The solid line shows the prescrip- 
tion bv lKewlev et al.1 (120011) . Objects with AGN-like spectra are 
close to the upper right while star-forming galaxies reside closer 
to the lower left corner as labelled. 



tions. Redshifts are derived from the I D spectra following the 
process described in Driv er et al. 



2.2. Emission line measurements 




0.1 0.2 
Redshift (z) 



Fig. 2: The redshift and absolute r-band magnitude distribution 
of galaxies in our 14 GAMA volume limited samples (see Table 
[TJ. White lines show the lower and upper apparent magnitude 
completeness of GAMA. 



We provide a brief summary of the spectroscopic process- 
ing, with full details provided by Hopkins et al. (in preparation). 
The blue and red parts of the AAOmega optical spectra are pro- 
cessed and r educed separately using the 2dfdr data re duction 
pipeline fsee lCroom etai1l2004t ISharp & Birchallll2010h before 
combination. The pipeline performs standard bias subtraction, 
flat-fielding, tram-line fitting and wavelength calibration using 
arcs. Sky subtraction is done using dedicated sky fibres before 
performing optimal ID extraction an d is refined using Principa l 
Component Analysis sky subtraction (Shar p & Parkinsonl2 010). 
This is followed by flux calibration using standard star observa- 



We use the Gas AND Abso rption Line Fitting algorithm 
(GANDALF. ISarzi et al.ll2006h . an extension of the Penalized 
Pixel cross-correlation Fitting (pPXF) algorithm of Cappellari & 
Emsellem d2004h . to model the stellar absorption and emission 
line spectral components of our flux calibrated spectra. We mea- 
sure the corrected flux and equivalent width of the strong emis- 
sion lines from the GANDALF output e mission line spectra. Our 
stellar population templates are from iMaraston & S tromback 
(2011) based on the c ode of Marastonj d2005h. which uses the 
MILES stellar library dSanch ez-Blaza uez et afll2006h with the 
revised spectral resolution Beifi ori et alJ d20i*ll) . Our emission 
line measurements are effectively stellar-absorption corrected 
through the template fitting performed by GANDALF. We use 
the Balme r decrement to apply the dust obscuration curve cor- 
rection of ICardelli et all dl989h to the emission line measure- 
ments, as recommended bv lCalzettild2001h . 



2.3. Sample selection 

We construct several su bsamples of the GAMA data. First, we 
use the relationship from Kauffma nn et al.1 (12003) to discriminate 
galaxies with significant AGN contribution fro m star-forming 
galax ies using their position in the BPT diagram dBaldwin et al.l 
1981). Fig. Q] illustrates this selection. Galaxies whose Ha or 
H/3 lines are affected by overlapping strong sky line residuals 
are rejected. We nominally select galaxies for which the strong 
emission lines Ha, H/3 and [NII]/16583 have S/N > 3. We refer 
to this selection as our "fiducial" selection throughout. A discus- 
sion of possible selection effects and justification of the fiducial 
selection are presented in Section 14.11 We reject any galaxies 
at z < 0.061 to minimise fibre aperture issues as described in 
IKE08I scaled for GAMA fibres (i.e., 2" in diameter). We also 
reject objects with a Balmer decrement Ha/H/3 < 2.5 and ex- 
tinction E(B - V) > 10. We refer to this sample as the 'main 
sample' for our analysis. 
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Finally, we also select many volume limited samples in nar- 
row redshift ranges, with maximum and minimum absolute r- 
band Petrosian magnitude determined based on the survey ap- 
parent magnitude completeness limits (i.e., 17.8 < r pe tro < 19.4). 
This is done in order to have the broadest range in absolute mag- 
nitudes possible, yielding a maximal stellar mass range in each 
bin. The redshift ranges for each bin are given in Table Q] and 
shown in Fig. [2] 



3. Analysis 

3.1. Star formation rates 

We compute the SFR from the stellar absorption and extinction 
corrected Ha flux measured by G ANDALF on the flux cali- 
brated spectra using equation 2 of iKennicuttl JT998). Because 
they are measured from the flux calibrated spectra, our SFR 
measurements include an implicit correction for fibre aperture 
effects. The implicit assumption made is that star formation is 
distributed across the galaxv following the stellar light (Hopkins 
et al. l2003h. Th roughout this work, we assume a Salneter IMF 
dSalpeterlll955h . O ur SFR measurements compa re well with 
those published in iGunawardhana et alj d2011h and SDSS- 
measured values in iBrinchmann et aTl d2004h . Throughout, SFR 
is quoted in units of M yr~' . 



tio. For N202 > -1.2 we solve the following polynomial: 

N202 = 1106.8660 -532. 1545 1Z + 96.373260Z 2 
-7.8106123Z 3 + 0.23928247Z 4 , 



(1) 



where Z=log(0/H)+12. 

If N2Q2 < —1.2. we use the average of the Kobulnicky & 
Kewlev d2004l hereafter KK041 and lMcGaughl dl99"Tl hereafter 
M91) methods. In M91, oxygen abundances (Zm9i) for the lower 
R23 branch are derived using: 

Z M 9i = 12-4.944+0.767x+0.602x 2 -;yi(0.29+0.332x-0.331x 2 ),(2) 



where yi 



log([Om]/U4959,5007/[On]/t3727) and x 



log/?23. KK04 obtain oxygen abundances (Zkkm) by iterating 
the following two equations: 

logq = (32.81 - 1.153y|+Z(-3.396-0.025y 2 (3) 

+0.1444$) x (4.603 - 0.31 I9y 2 - 0.163y 2 + Z 

-0.48 + 0.0271y 2 + 0.02037$)"', 

where y 2 = log([Om]i5007/[On]A3727) and 

Z KK 04 = 9.40+4.65x-3.17x 2 -log4(0.272+0.547x-0.513x 2 ),(4) 

with q the ionization factor, until convergence is achieved. 
Hence, for N202 < —1.2, the final quoted oxygen abundance 
is simply 12 + \og(0/H) = (Z M s>i + Z KK04 )/2. 



3.2. Stellar masses 



Stellar mass measurements are described in lTavlor et al.l (1201 lh . 
They are determined using spectral energy distribution fitting of 
the m-, g-, r-, i- an d z-magnitudes. We use a library of ~80,000 
templates from the iBruzual & Chariot] (|2003h stellar population 
models. The dust curve of lCalzettil d2001|) is assumed during this 
process. The quoted stellar mass (in units of M Q ) is derived from 
marginalisation over the posterior probability distribution of all 
fitted parameters (e.g., age, metallicity, stellar mass, star forma- 
tion history, etc). An offset of 0.2 is applied to the log of the 
stellar mass values published in iTavlor et al.1 (1201 lh to accom- 
modate our assumption of a Salpeter IMF. 

3.3. Gas-phase metallicities 

As stated in Section [1] several methods may be used to de- 
rive gas- phase metallicities of star-forming galaxies (see re- 
views byjLopez-Sanchez & Esteban KOlOklLopez -Sanchez et al.1 
|2012|) . Although the direct method should be preferred, the faint- 
ness of the auroral lines (particularly, the [OIII]/l4363 emis- 
sion line) often prevents its use for relatively high metallicities 
(12+log(0/H)>8.4). Hence, techniques using the more readily 
measurable strong emission-lines have been developed to esti- 
mate galaxy gas-phase metallicity. In this work, we present a 
representative set of three of the popular strong emission line 
methods discussed in the literatu re. These include te chniques 
based on photoionization models ( M9 it IKD02T: lKK04h and em- 
pirical calibrations from samples of objects with known metal- 
licity dPP04h . 

Method 1 : Kewlev & Dopital {2002) 



First, we use the method d escribe d in lKD02l with the update de- 
scribed in Appendix A of IKE08L This can be summarised as a 
combination of the N202 = [Nn]/16583/[Oii]^3727 ratio and 
an average of several calibrations of the degenerate R23 line ra- 



Method 2: lPettini & Paael |2004) 

The s econd method used is that described in iPettini &~ Pagel 
d2004L hereafter PP04). Briefly, we use the 03N2 ratio, which 
is defined as follows: 



(5) 



Q3N2 = log f[O"^006/ H/ n 
s \ [NH]A6583/Ha/ 

In |PP04, the 03N2 ratio abundances are calibrated using HII 
regions yielding the following relationship: 



12 + log(0/#) 



PPM 



8.73- 0.32 X03N2. 



(6) 



Because metallicities obtained directly using individual HII re 



gions result in lower metallicities (e.g. KD02), we use the public 



SDSS line measurement dBrinchmann et al 



2004) and oxygen 



abundances (TQ4J) t o der ive the following empirical conversion 
between PP04 and |T04| abundances (see Lara-Lopez et al., in 
preparation): 

12 + hg(0/H) = 0.103 + 1.021 x (12 + log(O/HW 04 ) ■ ( 7 ) 

Equation Q "accommodates" the well known 0.2-0.4 dex off- 
set expected between metallicities derived using the direct and 
model based methods (see Section[T]i. 

Method 3: lKobulnickv & Kewlev! d2004|) 

The third abu ndance determinati on met hod presented in this 
work is that of KK04 as updated by |KE08l It is based on the pop- 
ular R23 ratio and hence yields representative results of studies 
that use this ratio. The lower R23 branch conversion is given by 
iterating Equations[3]and|4]until convergence. For N202 > —1.2, 
metallicities are derived by iterating Equation [3] and 

12 + log(0/H) KKM = 9.72 - 0.777jc - 0.95 lx 2 - 0.072x 3 (8) 

-0.81 lx 4 - log(4)(0.0737 - 0.0713x - 0.141x 2 + 0.0373x 3 



-0.058x 4 ) 



until convergence. 



4 



C. Foster et al.: The mass-metallicity relationship in GAM A 




Fig. 3: A summary of the influence of varying the apparent selection criteria such as the survey apparent magnitude limit (r) and 
S/N limit in our fiducial lines (H a, Hfi and [NII]/16583) for three representative metallicity calibrations as labelled: KE08 (top row), 
PP04 (middle row) and KK04 (lower row). For each row, red contours show the 1, 2 and 3 cr profiles for our fiducial selection 
criteria, while the black contours show how varying S/N ratio cut and survey apparent magnitude limit influence the position and 
shape of the MMR. The red and black solid lines are respective second order polynomial fits to the MMR data. The fitted position 
of the MMR varies slightly for all ratios, but these variations are small (i.e., typically < 0.05 dex). 



Uncertainties 



4. Results 



In all cases, uncertainties are computed using Monte-Carlo 
methods. Assuming that the uncertainties on the measured emis- 
sion line fluxes are normally distributed, we draw 1000 random 
measurements of the relevant emission lines for each object and 
recompute the metallicity each time. The quoted uncertainty for 
each calibration is the standard deviation of the relevant 1000 
measurements. 



For completeness, we also compute oxygen abundances us- 
ing the methods described in IKD02I for various available emis- 
sion lines and the R23 calibrations o f IM9 lLIZaritskv et al J (1 1 994h 
and T04 . We choose to present the above mentioned methods of 
lKE08LlPP04l and the R23 calibration of lKK04l as a representative 
selection of the various methods used in the literature. We find 
that these three methods summarize and represent well the range 
of results obtained from the variety of measurement methods. 



This section focuses on the sensitivity of the fitted MMR with 
a set of common apparent selection criteria and abundance cali- 
brations. We also look for evidence of redshift evolution and de- 
pendency on SFR. The three-dimensional distribution of galaxy 
stellar mass, metallicity and SFR in the joint GAMA and SDSS 
samples is discussed in detail in a companion paper (Lara-Lopez 
et al., in preparation). 



4.1. Mass-metallicity relation 



As mentioned in Section I2.3I we study the dependence of the 
MMR on selection criteria such as the signal-to-noise (S/N) in 
various lines and apparent magnitude completeness, as well as 
on the metallicity calibration used. For this, we compare the 
MMR obtained with our fiducial selection criteria with that ob- 
tained by increasing the S/N limit from 3 to 8, but keeping the 
r-band Petrosian magnitude limit constant. Next, we vary the r- 
band Petrosian magnitude limit between 19.4 and 18.3, while 
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Fig. 4 : GAMA MMR measured using the a variety of selection criteria taken mostly from the recent lit e rature, as lab eled, for the 
PP04 abundance estimate. The MMR obtained using the selection criteria of Brinchmann et al. I (l2004h. IToH IKE08L Lara-Lopez 
et al. (|20104 LL10),|Kob ulnickv & Ke wlevl d2004 KK04) andl Mannucci et all d2010L M10) (respective black 1. 2 and 3cr contours) 
are compared to that obtained using the selection criteria used in this work (red 1, 2 and 3<x contours). The MMR measured varies 
significantly (> 0.05 dex) if one selects on the [OIII]/15007 line. 




1.5 9.0 9.5 10.0 10.5 11.0 11.5 12.0 8.5 9.0 9.5 10.0 10.5 11.0 11.5 12.0 8.5 9.0 9.5 10.0 10.5 11.0 11.5 12.0 

log M, 



Fig. 5: GAMA MMR measured using galaxies selected with metallicity uncertainties less than 0.1 dex (left), 0.05 dex (middle) and 
0.01 dex (right). Contours are similar to those of Fig. [3] and [4] 



keeping the S/N limit at 3. In all cases, we fit a secoi d or- Lopez etal 



der polynomial to the distribution obtained, following tit 
proach of other recent work. In what follows, we consider ver- 
tical (metallicity) variations of order 0.05 dex to be significant, 
as they could account for at least half of the reported variations 
of the MMR with redshift for z < 0.4 (see for example iLara- 



contours or 



2009b) and specific SFR (e.g.. lEllisonetal.l2008l) . In 



e ap- other words , we are looking for vertical differences of either the 



fitted MMR with variations of the selection. Fig. [3] 
summarizes the results of this exercise. Varying these selection 
criteria for Tie MMR does not significantly alter the fitted M MR 
for a given! calibration, however as emphasised by KE08, the 
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Fig. 6: The GAMA mass-metallicity (upper left), mass-SFR(upper right) and metallicity-SFR (lower-middle) relationships in our 
main sample. Black and coloured contours show 2<x contours for the main and v olume l imited samples wi th med ian redshifts as 
labelled, respectively. Red and black solid lines show the fitted SDSS MMR from lKE08l using the updated |KD02| conversion and 
our fit to the GAMA data, respectively. Open circles in the lower panel show the expected mean metallicity and SFR as inferred 
from the mean mass using Equations [9] and [10] for the lowest, intermediate and highest redshift bins. The actual mean metallicity 
and SFR for these samples are shown with crosses (see Section l4~3l l. 



shape and location of the MMR is significantly altered between 
the three calibrations shown. 

Fig. [4] shows that the MMR can be significantly different 
depending on the choice of selection criteria, for a variety of 
criteria that have been used in recent work. In general, we find 
that selecting on the [OIII]/15007 line most influences the posi- 
tion and shape of the measured MMR. This is a consequence of 
higher metallicities being associated with relatively fainter oxy- 
gen lines, which are more likely to show lower S/N. We note 
that the effect is more significant with the higher S/N threshold 
applied by KE08 and LL10, compa red to the lower threshold 
applied bv lBrinchmann et all d2004l) . although even this cut in- 
troduces a measurable effect. We therefore recommend caution, 
and rigorously consistent sample selection if cuts are being made 
based on the S/N ratio in these lines, when studying small varia- 
tions of the MMR. 

An alternative way to identify galaxies with reliable metal- 
licity measurements is to select based on the size of the metallic- 
ity uncertainty itself. The effects of such a selection are shown 



in Fig. [5] where we vary the selection on metallicity uncertainty 
from 0.1 to 0.01 dex. It can be seen that for reasonably mod- 
erate criteria, the MMR and contours are essentially unaltered, 
however applying an extreme cut does alter the shape and po- 
sition of the MMR. We cannot reject the possibility that these 
could be due to genuine evolutionary effects since galaxies at 
lower redshift would inevitably have lower errors. More impor- 
tantly, we insist that selecting based on overall metallicity uncer- 
tainties can only be applied for calibrations that are monotonic 
and smooth conversions between observ ables and metallicity. 
Piecewise calibrations such as the IKD02I and R23-based meth- 
ods have discontinuous conversion functions that yield large er- 
rors around the transition, hence a selection on error would dis- 
criminate against galaxies of those or similar metallicities (typi- 
cally 12 + \ogiPIH) ~ 8.4 dex). 

Fig. [6] shows the MMR, the mass-SFR relation and the 
metallicity-SFR relation for GAMA. Each panel shows the re- 
lation for our main sample in greyscale, with the distribution 
for each of the 14 volume limited redshift samples overlayed as 
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coloured contours. We emphasize that even though our volume 
limited samples are complete in a luminosity sense (absolute r- 
band magnitude), this does not correspond to a strict stellar mass 
limit. While each volume limited sample is itself complete, we 
observe the expected Malmquist-type bias as we consider sam- 
ples selected at different redshift. Galaxies of low (high) stellar 
masses are preferentially observed in the lowest (highest) red- 
shift samples. As a result, in each redshift bin, the range of stellar 
masses probed is limited. This complicates the study of possible 
evolution within GAMA. This limitation is addressed by com- 
bining the SDSS and GAMA samples by Lara-Lopez et al. (in 
preparation). 

Bearing all previously mentioned caveats in mind, we ten- 
tatively find that the MMR in GAMA is mostly consistent with 
that measured in SDSS. It is important to emphasize that even 
though GAMA metallicities are scaled to the same calibration as 
that used in the presented SDSS MMR they are not strictly mea- 
sured using the same selection criteria and metallicity diagnos- 
tic. Hence, for the reasons already outlined above, its exact po- 
sition and shape is not directly comparable to the GAMA MMR 
presented here. Nevertheless, the similarity between the SDSS 
and GAMA MMRs is encouraging to see given the inherently 
different apparent magnitude ranges probed by the two surveys. 
Indeed, although GAMA covers a smaller volume than SDSS, it 
is probing to fainter magnitudes, and a galaxy of a given mass 
will on average be detected at higher redshift in GAMA than in 
SDSS. This consideration suggests that any evolutionary effects 
over this redshift range are likely to be small. A companion pa- 
per (Lara-Lopez et al., in preparation) compares the SDSS and 
GAMA surveys directly and self-consistently. 

As is common, we fit a second order polynomial to the MMR 
in GAMA and obtain 

log(0///) + 12 = (-2.2 + 0.3) + (1.95 + 0.07) x log(M«) (9) 

-(0.085 ± 0.003) x log(M t ) 2 

with an rms s catter of cr = 0.12 dex, similar to that foun d 
for SDSS (e.g.. lTremonti et al.ll2004t iKewlev & Ellisonl l2008). 
Equation|9]is valid for 8.7 < log(M «(M )) < 12.5 for our fidu- 
cial selection criteria and the |PP04| metallicity calibration. We 
emphasize here that the MMR fit in this fashion arises from 
the contribution of galaxies at all redshifts probed in GAMA 
(0.061 < z < 0.4). As a result, galaxies at low (high) redshift 
generally contribute to the low (high) mass end of the MMR, and 
hence possible evolutionary processes are likely to be masked. 

4.2. Mass-SFR relation 

The upper right panel of Fig. [6] shows the relationship between 
stellar mass and SFR in GAMA for our main sample as well as 
the 14 volume limited redshift samples. While the relationship is 
difficult to detect in the individual volume limited samples due 
to the small stellar mass range probed in each one, we find a pos- 
itive correlation between log(M*) and SFR in the main sample, 
spanning the full redshift range probed. We use a bisector fit for 
galaxies in the main sample with reliable stellar mass and SFR 
measurements after applying a 2<x clipping. We obtain: 

log(SFR) = (1.304 + 0.008) x (log(M,)) - (12.98 + 0.08) (10) 

with a normal rms scatter of cr ~ 0.5 dex. As with the MMR, 
we find that all the volume-limited samples are consistent with 
the global fit, and the overall trend is consistent with that found 
in the SDSS, although the global fit is still affected by the 
Malmquist bias. 



4.3. Metallicity-SFR relation 

The lower panel in Fig.[6]shows the relationship between metal- 
licity and SFR for galaxies in GAMA. We find only a weak re- 
lation between the two in only some of our volume-limited red- 
shift samples. The main sample itself doesn't exhibit any corre- 
lation. In order to better understand the effect of the MMR and 
mass-SFR relationships measured in Sections [4. ll and l4~2l on the 
metallicity-SFR distribution, we perform the following exercise. 
We compute the mean stellar mass in each volume-limited sam- 
ple and use Equations l9l and [TOl to infer a corresponding metal- 
licity and SFR, respectively. These are shown for the lowest, in- 
termediate and highest redshift bins as open circles in the lower 
panel of Fig. [6] The actual mean metallicity and SFR for these 
samples are shown with crosses. While Equations l9l andfTOlpro- 
vide a consistent estimate for the position of the data at the in- 
termediate redshift samples, they seem to over- (under- )predict 
the metallicity at low (high) redshift by ~ 0.05 dex. This re- 
sult is likely to arise as a consequence of the fits to the whole 
sample being dominated by the intermediate redshift samples 
where the data density is greatest. The implication, though, is 
that at the higher-redshift end the observed metallicity is slightly 
lower than would be consistent with the extrapolation from the 
intermediate redshift objects. Conversely, at the lower-redshift 
end, the observed metallicity (and SFR) is slightly higher (~ 0.2 
dex). This may be suggestive of evolutionary effects, even over 
this redshift range (see also Lara-L0'pez et al., in preparation). 

5. Summary and Conclusions 

We have presented the dependence of the mass-metallicity re- 
lationship on a variety of selection criteria, metallicity calibra- 
tions and SFR. As in previous studies, we see the strong mass- 
metallicity relationship for the GAMA samples. The shape and 
position of the mass-metallicity relationship can vary signif- 
icantly depending on the metallicity calibration and selection 
used. In particular, we find that selecting on the [OIII]/15007 line 
most affects the position and shape of the mass-metallicity re- 
lationship, while varying survey depth does not. Alternatively, 
for monotonic metalllicity calibrations such as that of Pettini & 
Pagel 120041) . one can also select on the metallicity uncertainties 
itself without biasing the sample. Given these, caution must be 
used when comparing the mass-metallicity relationship obtained 
using different metallicity calibrations and/or sample selection 
criteria. 

The mass-metallicity relationship in GAMA is in reasonable 
agreement with that found in SDSS despite the very different 
stellar mass range probed at any given redshift. Notwithstanding 
the cautions illustrated in this work, this result implies that evo- 
lutionary effects out to z < 0.35 are likely to be small. Such small 
effects reinforce the importance of caution in selecting samples 
with which to estimate their measurement. 

Using volume-limited redshift samples, we find tentative ev- 
idence that galaxies at lower redshift may indeed have measur- 
ably higher metallicities than those at higher redshift, even over 
the redshift range probed within GAMA. This suggests that there 
could be detectable evolution of the MMR on the order of a few 
tenths of a dex in metallicity spanning 0.06 < z < 0.35. 

In a companion paper (Lara-Lopez et al., in preparation), we 
investigate the joint SDSS + GAMA mass-metallicity relation- 
ship, its evolution and the joint 3D distribution in metallicity- 
log(M»)-SFR space self-consistently and in detail. 

Finally, while the use of fibre-fed spectrographs has enabled 
large redshift surveys up to now, all are essentially limited by 
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aperture biases. Indeed, it is imperative to test the basic assump- 
tions generally used in similar work that 1) star formation fol- 
lows the distribution of stellar light and 2) models of single HII 
regions can be applied to spectra of entire galaxies. These as- 
sumptions can now be tested using available integral field unit 
spectrographs by observing galaxies at low redshift. Moreover, 
the Sydney-AAO Multi-obje ct Integral field sp ectrograph with 
its hexabundle technology ICroom et al.l d2012h will allow large 
spatially resolved spectroscopic galaxy surveys, enabling us to 
test these important assumptions and eliminating aperture bi- 
ases. 
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